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Aims. The Self-Similar Secondary Infall Model (SSIM) is modified to simulate a merger event. 

Methods. The model encompasses spherical versions of tidal stripping and dynamical friction that agrees with the Syer & White 
merger paradigm's behaviour. 

Results. The SSIM shows robustness in absorbing even comparable mass perturbations and returning to its original state. 
Conclusions. It suggests the approach to be invertible and allows to consider accretion as smooth mass inflow merging and mergers 
as intermittent mass inflow accretion. 
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1. Introduction 

Structure formation in the Cold Dark Matter (CDM, or 
more simply DM) paradigm is dominated by the hierarchi- 
cal pictu re of repeate d merg ers. This picture was empha- 
sised by ISver & White! dl998 ), explaining the dynamical for- 
mation of halo density profile with a feedback mechanism 
provided by repeated mergers. Whereas it is now believed 
that isotro pisation of the velocity d i spersi o n (angular momen- 
tum; see iLe Delliou & Henrikserj l2003h iBarnes et all 120051; 
MacMillan et al. 2006) via the radial-orbit insta bility (also 
viewe d as adiabatic variability of self-similarity, IHenriksenl 
120071) is responsible for the density profile formation, their pic- 
ture remains a widely accepted description of the merger diges- 
tion mechanism. Despite its simple spherical symmetry and ap- 
parent lack of compliance with the merger paradigm, some stud- 
ies have shown that the Secondary Infall Model (SIM) is a viable 
model to predict the structure and density profil e evolutions of 
DM haloes as compared to N-bo dy simulations ( Ascasiba r et alj 
l2007t ISaTvador-Sole et alJl2007h . 

This letter proposes to understand this paradox by examining 
the merger paradigm within the SIM and studying how merger 
events impact on the relaxation and structure of a CDM halo. 

Th e SIM stems from the seminal work of iGunn &Gottl 
dl972l), and the SS I M (S elf-si milar SIM) started when 



Fill more & Goldreichl d 19841) and iBertschingeil d 19841) inde 
pendently found self-similar solutions to the SIM. It was later 
shown that those solutions can be reached from non-self- 
similar initial condition s (e.g. in iHoffman & Shahaml 1 1985; 
White & Zaritskvl Il992t iRvdenl 11993) ; I Henriksen & Widrow l 



1995 



1999; 



119971; lAvila-Reeseet alJ Il999t IHenriksen & Widrov 
Idel Popolo et al J l2Q00b IHenriksen & Le Delliou! 12002; 
Le Delliou & Henri ksen 2003) an d a systematic approach to th e 
SSIM was used in IHenriksen & Widrowl dl995l [l997l 1 19991): 

d2002l) : ILe Delliou & Henrikser] 
the Carter-Henriksen formalism 



IHenriksen & Le Dellioul 
(2003^ derived from 



IWilliams et all l2004t iLu et alj 120061). b ut no fundamental 
attempt was made before Le DelliouT i 2002) to confront the SIM 
with the merger paradigm. 

The following section (Sec. |2| will describe how and why 
the SSIM can be extended to model a merger event. Then Sec. 
[3] will discuss how the symmetry of the SSIM still allows for a 
form of tidal stripping and dynamical friction, before presenting 
the consequences of such a merger in the SSIM in Sec.|H and to 
make some concluding remarks in Sec. [5] 



2. Merger in an Infall 

Modelling a merger event in a spherical geometry may appear 
contradictory but it is possible to a certain extent. To understand 
this it is important to realise the following: a ver y small amount 
of sub structures are seen in N-body simulations; Dieman d et al.l 
(2007) find that only 5 .3% of the total mass fraction of haloes 
lie in subhaloes. In the ISver& White! d 19981) picture, incoming 
satellite haloes merge with their parent, fall in the centre and 
contribute to the density profile and to the parent's relaxation and 
virialisation. However, in simulations, subobjects swing back 
and forth several times in their parents before being digested. 
That process can be modelled in a simpler way: on average, 
spherical symmetry is not bad dAscasibar et al.l 120071) as it re- 
produces the correct time scales and density profiles. Shell codes 
are much simpler than N-body codes and therefore provide with 
robust tests of certain aspects of their results. Other simplifying 
approaches have been used to u nderstand halo formation, such 
as phase-space coarse graining dLe Delliou & Henriksenl 120031: 
IHenriksenl 120041. 12 006) or in the one dimensional slab model 
used in Binney (2004), where it was shown to explain the for- 
mation of cosmic web sheets through the interplay of phase 
mixi ng and violent relaxation, a lso present in spherical mod- 
els. |Henriksen ^Widro\3 d!99 9) have shown that relaxation is 
moderately violent (in their figure 9) and induced by a phase 
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Table 1. Density, mass and mass perturbation ratios defining the 
satellite initial OD for the mergers in the SSIM. The mass per- 
turbation measures the perturbation of the OD compared to the 
background halo region it spans, just before entering the core. 
First column gives parent initial power law seed and panel order 
in reference to figures [TJ|2] and [3] 



In this pap er the SSIM is implemented in a shell code 
(see details in iLe Dellioul I2002L and references therein) with 
fully dynamical Lagrangia n treatment of infall using the CH 
(Carter & Henrikserj 1 19911) self-similar variables that reveals 
when the system reaches naturally a self-similar regime. A 
halo is modelled from a radial power law perturbation dp/p °c 
r~ e on an Einstein-de Sitter homogeneous background, that is 
evolve d to reach its quasi-stationa ry self-similar regime in its 
coreQ dHenriksen & Widrowlll999l). The SIM is known to es - 
tablish a self-similar infall phase (Henriksen & Widr owiri997l) . 
whic h then leads to a semi-universal power law d ensity pro- 
file dFillmore & Goldreich|[l984t lBertschingerlll984l) : for initial 
power index e < 2, the isothermal sphere (p oc r^' with p = 2) 
is the semi-universal attractor, whereas with e > 2, there is a 
continuum of attractors with p = 3e/(l + e). Positive overden- 
sity and the requirement of a finite initial core mass in the centre 
limit the range to < e < 3. The cores explored here were cho- 
sen, as presented in Table [TJ according to their SSIM behaviour 
defined by their initial power index: typical shallow (e = 3/2) 
and steep (e = 5/2) profiles, with the addition of an extreme 
steep case (e = 2.9) to test the behaviour of a highly concen- 
trated parent halo. The steep and shallow denominations refer to 
the comparison relative to the isothermal sphere. 

In this geometry, an overdensity (hereafter OD, or satellite), 
representing a spherically averaged satellite halo, is a region of 
overdense shells close to the edge of the core, the parent halo 
(hereafter core, or parent). 

The OD is evolved dynamically from an initial gauBian den- 
sity profile added on top of the background halo profile over a 
finite region. That evolution runs long enough to observe the sig- 
nature of the OD's own stationary regime in phase space. This is 
manifested in the mixing of its Liouville sheet during the OD's 
dynamical mass accretion of halo shells from its environment. 
The OD's definition as a set of particles (shells) is frozen when 
the core swallows it. 

At that point are recorded the ratios of OD-over-core masses, 
M ra ,i , of their densities, Z) ra/J0 , and the measure of the pertur- 
bation provided by the OD on its background surroundings, in 
mass, Mod/Mbg- For each case, three different satellites were 
chosen, trying to obtain various types of mass and density ratios 
between satellites and parents. 

Since they were allowed to accrete mass dynamically from 
their environment, ODs were laid close to the edge of the core 
to maintain some control over the final frozen mass and density 
ratios. Some configurations of those ratios were too difficult to 
obtain: in the shallow case, with high M ra , IO , lower values for 
D m tio were prevented by the high density background the OD 
accretes from, while for the steep cases, also with high M raflo , 
higher D mlio couldn't be obtained because of their cores' poor 
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Figure 1. Shallow case: Virial ratio and phase space diagrams 
at the end of the self-similar phase for three sets of ODs in the 
e — | = 1.5 case. Overdensity shells are emphasised in green 
in phase space. Digestion time is defined from OD core entry 
(pre-spike, see Sec. HI to virial recovery (measured on left pan- 
els). Zoomed encapsulation of those spikes in middle and lower 
left panels show their measure in troughless cases. T, X and Y 
are respectively the self-similar time, radius and radial velocity, 
which units are set by G — M(oo) = 1 dHenriksen & Widrowl 
119971) . 



density backgrounds which tended to spread the ODs (see Sec. 
Els tidal effect). 

The ratios indicated are measured at the time of core entry. 
The explored values are presented in Table Q] 

It is crucial to point out that the numerical implementation 
of the SSIM entails a shell code where finite size shells model 
the continuous system. That will play a role in the discussion of 
the results. 

3. Merger paradigm and SSIM 

ISver& White! (Il 998b have attempted to define the singularity of 
mergers in an effort, at the time, to explain the universality of 
the de nsity profile found in N-body simulation bv lNavarro et al.l 
( 1996, hereafter NFW): their key feature is the feedback mecha- 
nism between dynamical friction from the parent halo and tidal 
stripping of the satellite. Even though this is not anymore consid- 
ered to hold the key to the formation of the density profile, their 
merger digestion mechanisms is still widely accepted to describe 
the behaviour of satellites. I argue that both mechanisms can be 
modelled within the SSIM despite its one-dimensional nature. 

Tidal acceleration on an infinitesimal shell of mass dm — 
4npr 2 dr - located at radius r, containing the system mass M 
and with thickness dr - can be defined as the differential gravity 
between its boundaries. Defining the cumulative average density 
profile 

M(r) 



"47rr 3 /3' 



(1) 



the inward oriented elementary tidal acceleration reads, to lead- 
ing order, 



(2) 



dT =AnGdr\p--(p) r 



It is thus clear that regions of peak density below the cumulative 

/ _ - 2 / _\ \ :n : j* j_i i 
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Table 2. Digestion and dynamical times and strength parameter 
of the OD for the mergers in the SSIM. Again, first column gives 
parent initial power law seed and panel order in figures Q~| [2] and 

m 



Dynamical friction classically is defined as the creation of 
a wake by a moving mass in a gravitating medium which back 
reaction entails a net drag force upon the moving mass. In the 
SSIM, a massive shell is crossing the core's shell in its travelling 
inwards or outwards. We will see that a radial drag force, with 
the correct orientation, is experienced as a result of this motion 
in the spherically averaged model. 

This crossing of shells by the OD's results in shells just out- 
side of it feeling more or less mass pulling inwards, depending 
on the direction of the motion of the massive OD shells. That 
leads to a differential tightening or spreading of core's shell be- 
hind the moving mass, in the fashion of a wake. However in 
spherical symmetry, an outer wake does not contribute to the 
pull on the OD. Nevertheless, its mass corresponds to shells that 
defected from inside because of OD motion, and their effect can 
be seen in the dynamics (see AppendixlAli. 

In a similar fashion, the dynamical effect on the OD from its 
motion can be described in terms of a drag force: the crossing of 
core shells by massive OD shell lead to a decrease, or increase, 
of the resulting inner mass of the moving OD, depending on the 
direction of motion. Thus, with inner mass goes the inner pull, 
which can be interpreted a dragging force that adds to the total 
force, that should be experienced in the opposite direction of the 
motion. 

Therefore, the SSIM with an outer overdensity can be inter- 
preted to model the main features of the merger paradigm. 

4. Digestions 

Indeed, it is possible to keep track, in the Lagrangian shell 
model, of the defined satellite's (OD's) components once 
they have been absorbed by the parent (core). The core can 
be considered isolated at the end of the accretion phase 
dHenriksen & Widrowl l 19971) . The phase space configurations of 
simulated merged haloes are displayed on Figs. [1] [2] and [3js 
right panels, distinguishing between the core and OD's accreted 
shells. This reveals how the different ODs, in their various (shal- 
low or steep) environments, either retain some degree of coher- 
ence after being ingested by the core or have been digested and 
scattered over the core's phase space. 

The left panels of Figs.[TJ[2] and[3]examine the Virial ratios 
of the corresponding cores, and show a remarkable robustness 
in the SSIM: the quasi-stable self-similar phased is shown to be 
either marginally or strongly disturbed by the OD absorption, 
but to return to the original undisturbed level of the parent af- 
ter a digestion time T digestion, provided a mass flow still fuels the 
self-similar equilibrium. Digestion is manifested by a more or 
less pronounced initial decrease (entry of extra mass in core in- 
creases W), followed by a spike (first crossing of centre gives 
moD high velocities, thus peaks K) and then, for stronger dis- 



Virial ratio of core E=2 5 Phase Space with OD 
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Figure 2. Steep case: Virial ratio and phase space diagrams at 
the end of the self-similar phase, for three sets of ODs in the 
e = I = 2.5 case, including an emphasis on digested overdensity 
shells in phases space and a measure of digestion time. Same 
units as in Fig.Q] 



lower Virial. iHenriksen & Widrowlll999l) . Its deepness depends 
primarily on M ra; , . Digestion time measurements are shown on 
Figs. [H El and [3]s left panels (double horizontal arrows), and 
are summarised in Table [2] There, they are compared with the 
OD's free fall dynamical time through the core, Td y „amicai, also 
indicated on the figures. 7 'dynamical is defined as the free fall time 
to the centre of a test shell across a constant density distribu- 
tion, equivalent to the core, in self-similar variables. From Table 
[2] without Fig. [TJs two lowest panels, where the definition of 
T digestion is problematic, the average (j 'digestion II 'dynamical) = 3.33, 
with a standard deviation of 0.77, can be computed. It shows 
the core digests the OD in 2 to 4 passages in the central relax- 
ation region of phase space. This is comparable to the number of 
distinguishable Lagrange-Liouville streams present in the core's 
outer phase space regions, as seen from Figs. Q] [2] and[3]s right 
panels. 

From the OD's point of view, the mergers display their ef- 
fects in phase spaces, represented on Figs. Q] IU and [3]s right 
panels, on which two features are crucial: the spread (or com- 
pactness) of the OD over the core at the end of the infall phase 
and the presence of some, or all, of its shells in the centre of 
the core's phase space. This refl ects the digestion mechanisms 
adopted bv lSver & White! d 19981) . Their proposal aimed at a dy- 
namical explanation of the NFW profile. Although this expla- 
nation is not anymore considered (see Sec. Q}, it is interesting 
to note that the presently discussed single merger model in the 
SSIM shows signs of inflections (central flattening and edge 
steepening) from its semi-universal, almost isothermal, density 
profile. However this is not the focus of this paper. 

The OD's compactness resists to tidal stripping while its fi- 
nal presence in the centre is driven by dynamical friction. The 
fate of a model satellite in the SSIM displays behaviour well in 
agreement wit h the merger digestion mechanisms proposed by 
ISver & Whitet in the SSIM a combination of density and mass 
ratios leads to emphasise each effect. High Z) rafio s seem to be 
the dominant factor for OD's compactness, while high M ra , I0 s 
promote the sinking of the OD to the centre of the core's phase 
space. 

All possible qualitative types of behaviour are present: if 
both ratios, M rano and D rat - W , are strong enough, the OD survives 

_i a • a ..l c 1 /rr: m _ 1 1101 _ 
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Figure 3. Extreme steep case: Virial ratio and phase space dia- 
grams at the end of the self-similar phase, for three sets of ODs 
in the e = 2.9 case, including an emphasis on digested overden- 
sity shells in phase space and a measure of digestion time. Same 
units as in Fig.[T] 



OD around, but not reaching, the centre of phase space (Fig.HJs 
lower right panel). Finally if both ratios are too low, the OD is 
scattered without reaching the centre of phase space (Figs.[T]and 
|2|s middle and[3]s lower right panels). 

A step further in this phenomenology would be to note that 
a combination of both ratios should be taken (M rot!0 .D rat!0 , see 
Table |2), for which a threshold can be defined for reaching the 
centre and another for compactness of the OD. However this 
classification seems to require an additional dependency with the 
steepness of the initial profile. Indeed the available data offer dif- 
ferent ranges for each initial profile case. The shallow case calls 
for higher values for the M rat:o .D mt j thresholds than the steep 
cases. This reflects the shallow case's wider spread of material, 
compared with the steep cases, that the OD has to cross on its 
journey towards the centre of phase space. 

As an illustration of our model, we can assume the Milky 
Way (hereafter MW) to have a shallow profile and use the cor- 
responding reliable digestion time model, that is with e = 1.5, 
M rat i = 0.751 and J 'digestion - 2.50. The corresponding satel- 
lite S would have a mass Ms — 44Mlmc compared to the Large 
Magellanic Cloud (hereafter LMC), which is huge. The model 
then yields a very short digestion time, also compared with the 
age of the oldest stars in the MW Tmw = 13.2Gyr, as 



T di 



■gestum 



^584Myr 



Tmw 
22.6' 



(3) 



Its dynamical time T dynamical — 234Myr indicates that at the end 
of digestion, this satellite's shells would be lined between the 
second incoming and second outgoing inner streams of the core 
and the model suggests it to then sink to the centre by the end of 
the MW formation as seen on Fig. [Tfs upper right panel. 



5. Discussion and conclusions 

The SSIM has proven its capacity to model a merger event. Its 
simplicity allows one to probe the dynamics of the merger and 
the most remarkable result of this work shows that the self- 
similar quasi-stable regime of quasi-Virial equilibrium is ex- 
tremely robust to perturbations that can be of comparable size 



the satellite's particles across the parents and towards its cen- 
tre a gree with the tidal s tripping and dynamical friction picture 
from lSver & White! d 1998b . provided some adaptation to the lan- 
guage of the SSIM's symmetry. Finally, and this is the claim of 
this paper, the numerical implementation of the model requiring 
discretisation, the rapid oscillations of the Virial ratio in the ac- 
cretion phase offer a novel interpretation in the light of the SSIM 
merger model: instead of a continuous stream of mass, the model 
presents a repeated bombardment of finite mass shells that can 
be understood as small overdensities; Fig.QJs zoomed two low- 
est right panels show a spike to manifest the weakest mergers 
digestion; thus the wiggles in the Virial ratio can be interpreted 
as manifestation of repeated mergers that are at this level indis- 
tinguishable from accretion. Therefore there is no fundamental 
difference between mergers and accretion, the latter being a se- 
ries of repeated merger with vanishing mass, while the latter is 
just intermittent accretion . This reconciles approaches such as 
S alva dor-Sole et a il (120071) where accretion was presented as a 
memory loss mechanism, eliminating the need to refer to merg- 
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Appendix A: Spherical model of dynamical friction 

A thin OD shell in the inward/outward direction, adding mass 
±dnioD, crossing shells at r creates a differential acceleration 
w.r.t. the state without OD which induces an infinitesimal dis- 
placement, thus a wake, 

G(df) 2 dmoD 



dr ■ 



2r 2 



(A.l) 



This wake of mass dm w = prdr induces on the OD an acceler- 
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opposite to the direction of motion. In addition, the amplitude 
of the drag force is shown proportional to dm DP, related to 

M ra ti .D ratio- 



